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FIELD OF THE INVENTION 

The present invention relates to a method of determining the thermal profile of a 
drilling fluid in a well. 

During drilling, the mud injected into the drill string of the well and flowing back 
5 through the corresponding annulus undergoes great temperature variations. The fluid 
can encounter temperatures that can range from 2°C for deep offshore wells to more 
than 180°C for very hot wells. Many mud properties, such as rheology or density, 
depend on the temperature. Calculation of the pressure losses during drilling can 
therefore be improved if an estimation of the temperature profile in the well is known. It 
10 is therefore important to be able to predict the temperature profile in the flowing mud 
from well data and mud characteristics. 

Measurement of the thermal profile of the fluid in a well under drilling would 
require complete instrumentation of the well, i.e. installation of evenly spaced out 
detectors in the drill string and in the annulus, allowing temperature measurement at 
15 various depths. However, installing such a measuring system entails too many 
constraints ; only localized measurements picked up by devices mounted in the drill 
string allow to know certain temperature points on the path of the drilling fluid. 

BACKGROUND OF THE INVENTION 

In the face of this lack of data, analytic models based on heat transfer equations 
20 have been developed to evaluate the thermal profiles of the fluid along the well under 
drilling. Some of these analytic models are implemented in softwares and allow to 
provide an estimation of thermal profiles from a certain number of data more or less 



2 



difficult to obtain. Thus, knowing the characteristics of the site and of the drilling 
equipment, by giving a value of the ten^erature of the fluid at the well inlet, these 
softwares can predict the temperature profile of the drilling fluid. 

However, a comparison between the results obtained with analytic methods and the 
5 measurements obtained iti the field shows that there can be great differences. 
Furthermore, the complexity of softwares using numerical calculation methods makes 
real-time implementation thereof difficult. 

On the other hand, a study of the bibliography on thermal models shows a similarity 
in the form of temperature profiles in most cases, which turns on three points : inlet 
10 temperature, outlet temperatxire and bottomhole temperature. 

The aim of this study is thus to propose a method allowing real-time determination 
of a thermal profile in the mud from three measuring points available in the field, i.e. 
the injection temperature, the outlet temperature and the bottomhole temperature 
measured by a detector mounted on the string. The form of the profile between these 
15 three points is represented by a type curve representative of the thermal profiles in a 
well under drilling, estimated from physical considerations on thermal transfers in the 
weU. 

SUMMARY OF THE INVENTION 

The method of determining the thermal profile of a drilling fluid circulating in a 
20 well under drilling according to the invention is defined by the successive stages as 
follows : 
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a) determining a general expression 61 for the thermal profile of the fluid within the 
drill string in the well and a general expression 92 for a thermal profile of the fluid in 
the corresponding armulus, using the heat propagation equation that takes into account a 
thermal profile of the medium surrounding the well, 

5 b) measuring the temperature of the fluid at the well inlet, Tl, at the well bottom, T2, 
and at the well outlet, T3, 

c) laying down that expressions 61 and 92 meet the temperature boundary conditions 
Tl,T2andT3, 

d) drawing the thermal profiJeofthe drilling fluid as a fimction of the depth. 

10 In order to obtain, in real time, a temperature profile with the method presented 

above, stages b), c) and d) can be repeated. 

According to the method of the invention, in stage a), general expressions 61 and 62 
can comprise unknown constants, and in stage c), it can be laid down that expressions 
61 and 92 meet the temperature boundary conditions Tl, T2 and T3 by determining said 
15 unknown constants. 

In order to determine a general expression 61 for the thermal profile of the fluid 
within the drill string in the well and a general expression 62 for a thermal profile of the 
fluid in the corresponding annulus, it is possible, according to the method of the 
invention, to use in stage a) the heat propagation equation that takes into account at least 
20 the thermal equation of the medium surroxmding the well, the flow rate of the fluid and 
the balance of the thermal exchanges undergone by the fluid, said thermal exchanges 
comprising at least the exchanges between the ascending and descending drilling fluid, 



4 



and/or to use the equation of heat propagation in a homogeneous medium on a cylinder 
of infinite height centered on the well, said cylinder comprisiag the drill string that 
guides the descending fluid and the annulus arouad said drill string, which guides the 
ascending fluid. 

5 According to the method of the invention, general ejqjressions 81 and 92 obtained 

in stage a) can be split up into several independent equations and, in stage c), it can 
furthermore be laid down that the profiles and the derivatives of the thermal profiles of 
the fluid within the drill string and in the corresponding annulus are continuous. 

The method according to the invention can notably be used to calculate the pressure 
10 drops of the drilling fluid circulating in a well under drilling, or in another application, 
to determine the zones of hydrate formation in the fluid during drilling. 

In relation to the methods for determining the thermal profile of a drilling fluid in a 
well according to the prior art, the present invention notably affords the following 
advantages : 

15 - the temperature profile obtained is more accurate because it is determined from 
three drilling fluid temperature measurement points while keeping an analytic 
expression of the thermal profile between the measuring points which is physically 
justified, 

- the temperature measurements being performed all the time, the method allows to 
20 obtain the temperature profile in real time and to observe the evolution thereof with 
time. 
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BRIEF DESCRIPTION OF THE DRAWINGS 

Other features and advantages of the present invention will be clear from reading 
the description hereafter of non-limitative examples, with reference to the 
accompanying drawings wherein : 

5 - Figure 1 diagrammatically shows the architecture of a well under drilling, 

- Figures 2, 3 and 4 show the form of the temperature profile of the drilling fluid in a 
vertical onshore well, 

- Figxire 5 shows the form of the temperature profile of the drilling fluid in a vertical 
offshore well, 

10 - Figure 6 shows the form of the temperature profile of the drilling fluid in a deviated 
offshore well, 

- Figure 7 shows the evolution as a function of time of the tenq>erature profile of the 
drilling fluid in a vertical offshore well. 

DETAILED DESCRIPTION 

15 It is possible to give an analytic expression for the thermal profile in the well and 

the amrnlus by using quite simple heat exchange considerations, i.e. the heat 
propagation equation. 

This model is based on the establishment of the heat balances in the well. 
According to a first approach, only the steady states are considered (the drilling mud 
20 flow is assumed to be stabilized for some time so that the temperatures no longer 
evolve). Certain hypotheses are necessary for calculation : the heat exchanges are 
measured in a plane perpendicular to the laminar flow of the mud, the various constants 
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are assumed to be independent of the ten^rature, and fcially the influence of the 
ten^erature of the medium surrounding the weU shows on an apriori selected useful 
diameter Rf 

It is then sufScient to use the heat propagation equation in a homogeneous medium 
5 on a cylinder of infmite height centered on the weU shown in Figure 1. In each well 
section, the heat loss equality is written by considering two temperature functions : 0,(z) 
within the drill string and 82(2) m the annulus. 
Let 

9f be the temperature of the formation, 
10 ;wthe thermal conductivity of the medium surrounding the well, 

>ta the thermal conductivity of the tubing (metal), 

Cp the heat-capacity rate of the drilling fluid, 

Rl the inside radius of the drill string, 

R2 the outside radius of the drill string, 
15 Rt the radius of the annulus, 

Rf the effective radius (for heat supply) around the well, 

D the flow rate of the drilling fluid, 

p the density of the drilling fluid. 

The heat balances per unit of depth are as follows : 
20 - Heat supphed by the medium surrounding the well to the fluid in the annulus : 
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- Heat carried from the fluid in the annulus to the fluid within the drill string : 

- Heat accumulated by the fluid in the drill string and in the annulus : 
Q, = -D.p.CpAGi 

Qa = D.p.CpA92. 

The heat balances lead to the following system : 
Qt = Q2 

Qa=Ql + Q2, 



These equations are solved by diagonalization and matrix inversion , and they lead 
10 to the folio wmg results : 



e^{z) = -K,(B + ry'-' - K,(B + r^Y^-' + 9, 
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with : 




2m, 



A = 



B = 



A + ^A^+AAB 
2 



A-^A^ + AAB 
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0f = a.z-i-6o being the thermal equation of the medium sxirrounding the well and a the 
thermal gradient. 

5 Ki and are the integration constants depending on the boundary conditions. 

It is thus possible, by means of some simplifying hypotheses, to obtain an analytic 
expression for the temperature profile of the drilling fluid in a well. If all the parameters 
are known, by giving the inlet temperature and by writing that temperatures 0i and 02 
are equal at the weU bottom, the profile is entirely determiaed. The most commonly 
10 used softwares use this type of predictive procedure. However, a study of the results of 
the models compared to field data shows how difficult it is to use these estimations in a 
predictive way. 

In the present invention, the system is based on the knowledge of three measuring 
points in the field : inlet temperature, outlet temperature and bottomhole temperature. In 
15 order to estimate the thermal profile in the well fi-om the three measurements consisting 
of the surface injection and outlet temperatures and the bottomhole temperature (inside 
or outside the drill string), the method according to the invention consists in coimecting 
the three measuring points by a general expression representative of the evolution of a 
thermal profile in a wellbore, as obtained according to the method described above. 
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We therefore take the eqiiations obtained by means of heat exchange calcvilations : 

According to the invention, these curve forms are adjusted to the three measuring 
points of the drilling fluid temperature at the inlet, Tl, at the well bottom, T2, and at the 
well outlet, T3. In order to use these three measuring points as boundary conditions, we 

5 choose to decouple the two equations (in the drill string and in the annulus) by using 
different integration constants while keeping the general expression. We obtain two 
general expressions of the temperature profile in the drill strmg, 61, and in the annulus, 
92, which have a physical significance but which comprise two degrees of freedom. 
Thus, ejq)ressions 01 and 02 can be adjusted by fixing said degrees of freedom in order 

10 to meet the temperature conditions Tl, T2 and T3. We therefore decide that the 
equations in the pipes and in the annulus have the form 2& follows : 

9, {z) = -K,Be'' ' - K^Be'^' + 0^-- 
B 

e^{2) = -K,(B + r,)e'''' - K^B + r^)e''^-' + 6^. 

We thus have four integration constants Ki, K2, K3 and K4 rather than two, which 
requires four boundary conditions to determine the temperature profile. These four 
boundary conditions are : measurements of the well inlet, bottomhole and outlet 
15 temperature, and an equality condition at the bottom between the temqperature in the 
driU string, 91, and the temperature in the annulus, 02. The profile can be adjusted at 
any time to pass through the measuring points : we thus have a real-time estimation of 
the thermal profile. Programming by means of a spreadsheet type software allows to 
readily obtain the real-time evolution of the profile. 
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Figiores 2, 3 and 4 respectively show the temperature profile of the drilling fluid in a 
vertical onshore well at a flow rate of 500 1/min, 1000 I/min and 2000 1/min. The 
analytic expression determined allows simple calculation of the temperature T in 
Celsius degrees of the fluid in the drill string (curve 91) and in the annulus (curve 92) as 

5 a fimction of the depth P in meter. The analytic expression depends on several 
parameters that can be fixed from the start. We use by default typical values of these 
parameters. To determine the terrperature profile of Figures 2, 3 and 4, the geothermal 
gradient a is assumed to be constant to correspond to the onshore situation of the well. 
The tercqperature profile is entirely determined by measuring the temperature, 20°C at 

10 the inlet, 35°C at the bottom and 24°C at the outlet of the well. 

The case of the vertical offshore well can be tackled by considering that the 
geothermal profile of the medium siirrounding the well is divided in two domains : let 
Om be the thermal profile of the sea and 0s the thermal profile of the ground. Thermal 
gradient a is assumed to be constant in each domain, but discontinuous from one 

15 domain to the other. Let am be the thermal gradient of the sea and as the thermal 
gradient of the ground. We then consider two series of equations (one for each domain) 
for each general expression in the pipes and in the annulus. We thus obtain four 
decoupled equations which represent the thermal profile of the drilling fluid in the well. 
Equation 011(z) corresponds to the temperature profile in the drill string in the sea, 

20 912(z) corresponds to the temperature profile in the drill string in the ground, 921(z) 
corresponds to the temperature profile in the annulus in the ground and 022(z) 
corresponds to the temperature profile in the annulus in the sea, 911 beii^ independent 
of 912 and 021 being independent of 922 : 



11 



6*12 (z) = -K^Be"'-' - K^Be'^ ' + -— 
B 

(z) = -K,iB + r{)e'' ' -K^(B + r^)e'^-' + 6, 
e^{z) = -K, (B + r,y'-^ -K,(B + r^)e^'-^ + 0^ 

This brings the number of integration constants to eight (Kj to Kg). The boundary 
conditions are then : measurements of the inlet, outlet and bottomhole temperature, 
equality condition at the bottom between the pipe temperature and the annulus 

5 temperati3re, to which we add the continuity of the thermal profiles in the drill string 
and in the annulus at the junction of the two domains and the continmty of the 
derivative of the thermal profiles in the drill string and in the annulus at the junction of 
the two domains. Similarly, it is then possible to obtain in real time a physically realistic 
thermal profile passing through the measuring points. Figure 5 shows the thermal 

10 temperature profile of a drilling fluid in an offshore well ftom the four equations 011, 
612, G21 and G22. The fluid circulates at 500 1/min and the temperatures measured are 
20°C at the inlet, 15°C at the outlet and 30°C at the well bottom. The thermal gradients 
are selected constant in each domain crossed by the well. 

Deviated wells represent the majority of the current wellbores. The physical 
15 problem is not fimdamentally different and it can be handled in the same way as 
offshore wellbores : the well just has to be divided into two domains, each domain 
being characterized by a different thermal gradient corresponding to the medium 
surroimding the well. In the case of a deviated well, the depth corresponds to the 
distance covered along the well trajectory. General expressions 01 and 02 representative 
20 of the thermal profile are each split up into two independent eqimtions. The vertical part 
is characterized by the thermal gradient a of the medium surrounding the well, the 
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deviated part is characterized by an eqiiation of the thermal profile of the medium 
surroxinding the well 9d = a.sin(<|)).z+0o, ^ being the angle of inclination. The same 
boundary conditions (measurements of the temperature at the inlet, at the outlet and at 
the bottom of the well, equality at the bottom between the pipe temperature and the 
5 annulus ten^erature, and continuity of the thermal profiles and of the derivative of the 
thermal profiles in the drill string and in the annulus at the junction of the two domains) 
then allow to solve the equations and to obtain the expression of the temqperature profile 
in the pipes and in the annulus. 

It is possible to combine the procedure applied for a vertical offshore well and the 
10 procedure applied for a deviated onshore well in order to determine the temperature 
profile in an offshore well when the direction of the hole is deviated in the ground. The 
domain is divided into three different domains : let 9m be the thermal profile of the 
vertical domain in the sea, 0s the thermal profile of the vertical domain in the ground 
and Gd the thermal profile of the deviated domain in the ground. Figure 6 shows the 
15 thermal profile in a deviated ofehore welL The fluid circulates at 500 I/min and the 
temperatures measured are 20°C at the inlet, 23°C at the bottom and 15°C at the outlet 
of the well. 

According to the same method as that used for the vertical offshore well or for the 
deviated onshore well, it is possible to determine the thermal profile of a vertical 
20 onshore weU whose formation thermal gradient changes as a fimction of the depth. The 
well is divided into domains characterized by a thermal equation of the medium 
surrounding the well. General expressions 61 and 92 representative of the thermal 
profile are then each split up into as many independent equations as there are different 
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domains. The same boundary conditions (measurements of the temperature at the inlet, 
at the outlet and at the bottom of the well, equality at the bottom between the pipe 
temperature and the annulus ten^rature, and continuity of the thermal profiles and of 
the derivative of the thermal profiles in the drill string and in the annulus at the junction 
5 of the two domains) then allow to solve the equations and to obtain the expression of 
the temperature profile in the pipes and in the annulus. 

By repeating the calculation allowing to obtain the expression of the temperature 
profile of the drilling fluid upon each new temperature measurement, we obtain a 
representation of the ten^erature profile evolving with time. Figure 7 shows the 

10 evolution of the temperature profile of the drilling fluid in an offshore well in the course 
of time. The graph in the up|>er part of Figure 7 shows the evolution as a fimction of 
time t in second of the flow rate parameter D in Vtnin of the drilling fluid, and of the 
temperature parameter T in °C of the drilling fluid at the inlet, Tl, at the bottom, T2, 
and at the outlet, T3, of the well. The three graphs in the lower part show the 

15 ten^rature profile at three different times and allow to observe the evolution thereof. 

Knowledge of the thermal profile of the drilling fluid at any time allows real-time 
calculation of the presswe drops in the well by taking into accoimt the thermal effects. 
This gives a better estimation of the bottomhole pressures and of the injection pressures 
for complex weUs. 

20 Another use of real-time determination of the thermal profile of the drilling fluid is 

hydrates formation prevention. Hydrates form under low temperature and high pressure 
conditions, conditions which are met notably in deep o:^hore weUs at the ground/sea 
interface. Knowledge of the temperature profile allows to determine the zones where the 
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temperature of the drilling fluid is below the minimmn value from which hydrates form, 
then to react accordingly, for example by raising the flow rate or by heating the fluid in 
order to prevent this formation of hydrates. 



